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Proper Motion Dispersions of Red Clump Giants in the Galactic 
Bulge: Observations and Model Comparisons 
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ABSTRACT 

Red clump giants in the Galactic bulge are approximate standard candles and hence they can 
be used as distance indicators. We compute the proper motion dispersions of RCG stars in 
the Galactic bulge using the proper motion catalogu e from the sec ond phase of the Optical 
Gravitational Microlensing Experiment (OGLE-II. ISumi et alJ20 04) for 45 fields. The proper 
motion dispersions are measured to a few per cent accuracy due to the large number of stars 
in the fields. The observational sample is comprised of 577736 stars. These observed data are 
compared to a state-of-the-art particle simulation of the Galactic bulge region. The predictions 
are in rough agreement with observations, but appear to be too anisotropic in the velocity 
ellipsoid. We note that there is significant field-to-field variation in the observed proper motion 
dispersions. This could either be a real feature, or due to some unknown systematic effect. 

Key words: gravitational lensing - Galaxy: bulge - Galaxy: centre - Galaxy: kinematics and 
dynamics - Galaxy: structure 



1 INTRODUCTION 



Many lines of evidence suggest th e presence of a bar at the Galac- 
tic ce ntre, such as infra red maps jDwek et alj 19951: Binnev et al 
1997) and star counts ( Stan ek et al. 1 1 19971 ; Nikolaev & Weinberg 



1997l : IUnavane & Gilm ore 1998), see I Gerhard ( 12003) for a review. 
However, the bar parameters are not well determined. For example, 
recent infra-red star counts collected by the Spitzer Space Tele- 
scope are best explained assuming a bar at a ~ 44° angle to the 
Sun-Galactic centre line ^Benjamin et al.|[2005h while most previ- 
ous studies prefer a bar at ~ 20°. In addition, there may be some 
fine features, such as a ring in the Galactic bu lge, that are not yet 
firmly established (Babusia ux & Gilmorel2005h . It is therefore cru- 
cial to obtain as many constraints as possible in order to better un- 
derstand the structure of the inner Galaxy. 

Many microle nsing groups monito r the Galact ic bulge, in- 
cludin g the EROsJXubourg et al]ll993l). MACHO fAlcock et al.1 
l200d). MQA teond et al.1 l200ll : fsumi et alj l2003tj) and OGLE 
1 Udalski et alf2 000) collaborations. In addition to discovering mi- 
crolensing events, these groups have also accumulated a huge 
amount of data about the stars in the Galactic bulge spanning sev- 
eral years to a decade and a half. 
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lEver & Woz niak (2001) first demonstrated that the data can 
be used to i nfer th e proper motions of stars, down to ~ masyr -1 . 
ISumietalJd2004l) obtained the proper motions for millions of stars 
in the OGLE-II database for a large area of the sky. In this pa- 
per, we focus on the red clump giants. These stars are bright and 
they are approximately standard candles, hence their magnitudes 
can be taken as a crude measure of their distances. As the OGLE-II 
proper motions are relative, in this paper we compute the proper 
motion dispersion s of bulge stars for all field data presented by 
ISumi et alj ( |2004|) . as they are independent of the unknown proper 
motion zero-points. These results could aid theoretical modelling 
efforts for the central regions of the Galaxy. 

The structure of the paper is as follows. In section [2] we 
describe the OGLE-II proper motion catalogue and compute the 
proper motion dispersions for bulge stars in 45 OGLE-II fields. In 
section [3] we describe the stellar-dynamical model of the Galaxy 
used in this work and detail how the model was used to generate 
proper motion dispersions. These model predictions are compared 
to the observational results in section|4]and in section|5]we discuss 
the implications of the results. 



2 OBSERVED PROPER MOTION DISPERSIONS 

The second phase of the OGLE experiment observed the Galactic 
Centre in 49 fields using the 1.3m Warsaw telescope at the Las 
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Figure 2. Extinction-corrected colour-magnitude diagram for stars in the 
OGLE-II field 1. The ellipse defines the selection criteria for RCG stars 
based on colour and magnitude, see text. Sample stars are also required to 
have proper motion errors si b < 1 masyr -1 and total proper motion 

< 10 masyr -1 . 



proper motion fti with measurement error is: 
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Maximising the likelihood ln(L) = ln(f|pi) for fi and a over all 
observations we find: 



d\nL 

djl 



= E 



(Mi - m) 
^ 2 + C 2 



= 



E 



and 



dlnL 



= E 



* 2 +c 2 



E 



(Mi ~ M) 2 
(* 2 +0 : 



= 



(2) 



(3) 



(4) 



which can be solved numerically to find a . 

The values of M and a obtained using the above maximum- 
likelihood a nalysis are virtually ident ical to those obtained via the 
equations in lSpaenhauer et al ] h992h . The errors on the observed 
proper motion dispersion values were determined from a bootstrap 
analysis using 500 samplings of the observed dataset. 



Campanas Observatory, Chile. Data were collected over an interval 
of almost four years, between 1997 and 2000. Each field is 0.24° x 
0.95° in size. Fig. Q] shows the position of the OGLE-II Galactic 
Bulge fields which returned data used in this paper. 



2.1 Red Clump Giants 

The red clump g iants are metal-rich horizontal branch stars 
dStanek et al.ll2000l and references therein). Theoretically, one ex- 
pects their magnitudes to ha ve (small) variations wi th metallicity, 
age and initial stellar mass dGirardi & Salarisll200lh . Empirically 
they appear to be reasonable sta ndard candles in the J-band wit h 
little dependence on metallicities dUdalskil2 000: Zhao et alj200lh . 
Below we describe the selection of RCG stars in more detail. 



2.2 OGLE-II proper motion data 

Bulge RCG stars are selected from the OGLE-II proper motion cat- 
alogue by applying a cut in magnitude and colour to all stars in each 
of the OGLE-II fields. We corrected fo r extinction and reddening 
using the maps presented bv lSumil f2004) for each field. Stars were 
selected which are located in an ellipse with centre (V — I)o = 1-0 
, Io = 14.6; and semi-major (magnitude) and semi-minor (colour) 
axes of 0.9 and 0.4 respectively, see Fig. [2] a similar selection cri- 
terion was used by Sumi (2004). Stars with errors in proper mo- 
tion greater than 1 masyr -1 in either the I or b directions were 
excluded. Stars with total proper motion greater than 10 masyr -1 
where similarly excluded, as these are likely to be nearby disk stars, 
see also section [372] Fields 44, 47-49 were not analysed due to the 
low number of RCG stars appearing in these fields. 

The proper motion dispersions for the longitude and latitude 
directions (o\ and crt,) were comp uted for each field via a maxi- 
mum likelihood analysis following Lu pton et alj dl987h . Assuming 
a Gaussian distribution of proper motions with mean jl and intrin- 
sic proper motion dispersion a, the probability of a single observed 



2.3 Extinction 

In order to ensure the correction for extinction and reddening above 
does not affect the kinematic measurements, <j\ and Ob were recom- 
puted for each OG LE-II field u s ing re ddening-independent mag- 
nitudes. Following IStanek et alj J 19971) we define the reddening- 
independent magnitude Iv-i- 



Jv-i = / - Al/ {Aw - Ai) (V - I) 



(5) 



where Ai and Ay are the extinctions in the / and V bands deter- 
mined bv lSumil fe004l) . The position of the red clump in the Jv-i, 
(V — I) CMD varies from field to field. The red clump stars were 
extracted by iteratively applying a selectio n ellipse computed from 
the moments of the data dRocha et alj 20 021 rather than centred on a 
fixed colour and magnitude. The selection ellipse was recomputed 
iteratively for each sample until convergence. The proper motion 
dispersions <ti and <Tb computed using RCG stars selected in this 
way are consistent with those determined using the original selec- 
tion criteria on corrected magnitudes and colours. 



2.4 Results 

Table Q] lists the observed proper motion dispersions along with 
errors for each of the 45 OGLE-II fields considered in this paper. 

Figures|3]and|4]show the proper motion dispersions o\ and <7b 
as a function of Galactic longitude and latitude. A typical value of 
ai or <7b of 3.0 masyr -1 corresponds to ~ 110 kms -1 , assum- 
ing a distance to the Galactic centre of 8 kpc. The proper motion 
dispersion profiles as a function of Galactic longitude shows some 
slight asymmetry about the Galactic centre. Thi s asymmetry ma; 



be related to the tri-ax i al Galactic bar structu re (Stanek et 



letry may 
t aljl997t 



iNishivama et alj2005l:lBabusiaux & Gilmorell2005l) . The most dis- 
crepant points in Fig. [3] correspond to the low-latitude fields num- 
bers 6 and 7 (see Fig. QJ. The varying field latitude accounts for 
some of the scatter in Fig. [3] however we note below in section |4~T| 
that there are significant variations in the observed proper motion 
dispersion between some pairs of adjacent fields. Owing to the the 
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Figure 1. The position of the 45 OGLE-II fields used in this analysis. The field used in Spaenhauer et al. 11992) is shown, located within OGLE-II field 45 
with (I, b) = (1.0245°, -3.9253°). 



lack of fields at positive Galactic latitude, any asymmetry about the 
Galactic centre in the proper motion dispersions as a function of 
Galactic latitude is not obvious, see Fig. [4] Field-to-field variations 
in the proper motion dispersions similarly contribute to the scat- 
ter seen in Fig. [4] along with the wide range of field longitudes, 
especially for fields with —4° < b < —3°. 

Table [2] lists the proper motion dispersions and cross- 
correlation term C*ib in the OGLE -II Baade's Window fie lds 45 
and 46 along with those found by iKoztowski et alj d2006t) using 
HST data in four BW fields. The two sets of proper motion dis- 
persions results are consistent at the ~ 2a level. It is important to 
note that the errors on the proper motion dispersions in Table[TJdo 
not include systematic errors. We also note that the selection cri- 
teria applied to stars in the HST data are very different to those 
for the ground-based data, in particular the magnitude limits ap- 
plied in each case. Th e bulge kinematics from the HST data of 
lKozIowskietalJ ( l2006h were determined for stars with magnitudes 
18.0 < /f814w < 21.5. The approximate reddening-independent 
magnitude range for the OGLE-II data was 12.5 < 7v-i < 14.6. 
The effects of blending are also very different in the two datasets. It 
is therefore very reassu ring that our resul t s are in general agreement 
with those obtained by Kozlows ki et al.1 j2006h using higher reso- 
lution data from the HST. For more comparisons between ground 
and HST RCG proper motion dispersions, see section|4] 

Figure [5] shows the cross-correlation term Cib as a function of 
Galactic co-ordinate. There is a clear sinusoidal structure in the Cib 
data as a function of Galactic longitude, with the degree of corre- 
lation between a\ and <7b changing most rapidly near I ~ 0°. The 
Cib data as a function of Galactic latitude may also show some ev- 
idence of structure. It is possible however, that this apparent struc- 



ture is due to the different number of fields at each latitude, rather 
than some real physical cause. 



3 GALACTIC MODEL 

The stellar-dynamical model used i n this work was produc ed us- 
ing the made-to-measure method JSver & Tremain el 1 19961) . The 
model is constrained to reproduce the density distribution con- 
structed from the d ust-corrected L-band COBE/DIRBE map of 
Sper gel et al.1 {T?96). An earlier dynamical model was built to 
match the total column density of the disk l|Bissantz & Gerhardl 
2002). This dynamical model matched the radial velocity and 
proper motion data in two fields (including Baade's window) quite 
well. No kinematic constraints were impose d during the construc- 
tion of the model. We refer the readers to iBissantz et alj ( 2004) 
for more detailed descriptions. The model used here is constructed 
as in that case with the further refinement that the vertical den- 
sity distribution is also included. This is necessary as the vertical 
kinematics (<7b) will also be compared with observations in this 
paper. However the density distribution near the mid-plane is con- 
siderably more uncertain, in part because of the dust extinction cor- 
rection. Thus the model used in this paper can only be considered 
illustrative, not final. Further efforts to model the vertical density 
distribution are currently under way and will be reported elsewhere 
(Debattista et al. 2007, in preparation). 

In Fig. [6] we present the mean motion of stars in the mid-plane 
of the Galaxy from this model. A bar position angle of 6 = 20° is 
shown here, as this is the orientation favoured both by optical depth 
measurements jEvans & Belokurovll2002l) and by the red clump gi- 
ant brightness distribution i Stanek et alj 19971) and was the angle 



4 Rattenbury et al. 



Table 1. Observed proper motion dispersions in the longitude and latitude directions, a\, tr b , and cross-correlation t erm Cp, for bul ge stars in 45 OGLE-II 
fields. High precision proper motion data for bulge stars were extracted from the OGLE-II proper motion catalogue <Sumi et al ] |20M) . N is the number of 
stars selected from each field. Fields 44, 47-49 were not analysed due to the low number of RCG stars appearing in these fields. 
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Table 2. Comparison between proper motion dis persions and cross-corr elation term Cib m two °f me OGLE-II fields (45 and 46) with proper motion 
dispersions computed from four nearby HST fields (Kozlowski et al. 2006). 
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1 Kozt owskietal. (2006) 2 This work. 
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Figure 3. Proper motion dispersion in the Galactic longitude and lati- 
tude (<7b) directions for 45 OGLE-II Galactic bulge fields as a function of 
field Galactic longitude. Open circles correspond to fields 6, 7, 14, 15 and 
43 which have relatively extreme galactic latitudes, see Fig.[T] 
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Figure 4. Proper motion dispersion in the Galactic longitude (<tj) and lati- 
tude (<Th) directions for 45 OGLE-II Galactic bulge fields as a function of 
field Galactic latitude. Open circles correspond to fields 6, 7, 14, 15 and 43 
which have relatively extreme galactic latitudes, see Fig.[T] 



used in deriving the model. Clearly one can see that the mean mo- 
tion follows elliptical paths a round the Galactic b ar. The analysis of 
OGLE-II proper motions bv lSumi et al. (2003b) is consistent with 
this streaming motion. 

3.1 Model stellar magnitudes 

The model has a four-fold symmetry, obtained by a rotation of 7r 
radians around the vertical axis and by positioning the Sun above 
or below the mid-plane. The kinematics of model particles falling 
within the solid angle of each OGLE-II field were combined to 
those from the three other equivalent lines-of-sight. This procedure 
allows an increase in the number of model particles used for the 
predictions of stellar kinematics. 

We assign magnitudes to stars in the Galactic model described 
above which appear in the same fields as that observed by the 
OGLE collaboration. Number counts as a function of /-band ap- 
parent magnitude, /, were used to compute the fraction of RCG 
stars in each of the OGLE-II fields. Figure [7] shows an example of 
the fitted number count function Nk (I) for one of the k = 1 . . . 49 
OGLE-II fields, wher e Nk (I) is of the form of a power-law and a 
Gaussian ( Sumi 2004J): 



N k (I) 



a fe 10^ + c fe exp 



2*1 



(6) 



where the constants a^, bk, Ck, I p> k,&k are determined for each of 
the k OGLE-II fields, see Tabled The fraction R k of RCG stars 
is evaluated as the ratio of the area under the Gaussian component 
of equation <(6j> to the area under the full expression. The integrals 
are taken over ±3ofc around the RCG peak in Nk(I) for each of 
the k OGLE-II fields. Fields 44 and 47-49 are not included as there 
are insufficient RCGs in the OGLE-II fields to fit equation (6). Fig- 
ure [7J shows that the model number count function fails to fit the 
observed number counts well for magnitudes / ~ 15.4. In order to 
convert stellar density t o a distribution of apparent magnitude, the 
relevant quantity is pr 3 l lBissantz & Gerharq|2002h . Depending on 
the line-of-sight, this quantity can give asymmetric magnitude dis- 
tributions through the bulge. Using the best-fitting analytic tri-axial 
density models for the bulge (Rattenbury et al. 2007, in prepara- 
tion), this asymmetry is observed and may explain the excess of 
stars in the number count histograms, compared to the best-fitting 
two-component fit of equation ^6). The inability of equation l|6} to 
model completely all features in the observed number counts in 
some cases leads to an additional uncertainty in the magnitude lo- 
cation of the fitted Gaussian peak. Computing the apparent magni- 
tude distribution as oc pr A also produces a small shift in the peak of 
the magnitude distribution. This shift is ~ +0.04 mag for 1 — 0°, 
b = 0°. The proper motion dispersions computed here are unlikely 
to be sensitive to these small offsets. 
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Table 3. Values of fitted parameters in equation (6) for all 45 OGLE-II fields used in this analysis. R is the ratio of observed RCG stars to the total number 
of stars in each field, evaluated over ±3<r around the RCG peak magnitude, I p , where a is the fitted Gaussian spread in equation {6}. The magnitudes of 
the model RCG stars are shifted by Am to correspond with the observed mean RCG magnitude in each field. The total number of model stars in each field 
assigned RCG magnitudes and colours is n rcg and the total number of model stars in each field is n a n. The corresponding total model weight values for each 
field are given by u> r cg and «j a n respectively. The large values of <r for fields 8-11 might be related to their position at large positive longitudes, and could 
indicate a structure such as the end of the bar, a ring or spiral arm. An analysis of the bar morphology based on these results is underway (Rattenbury et al. 
2007, in preparation). 
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Each star in the galactic model is assigned a RCG magnitude 
with probability Ru for each field. The apparent magnitude is com- 
puted using the model distance. Stars which are not assigned a RCG 
magnitude are assigned a magnitude using the power-law compo- 
nent of equation ^6), defined over the same limits used to compute 
Rk . Here we implicitly assume that the RCG stars trace the overall 
Galactic disk and bulge populations. 

The RCG luminosity function is approximated by a Gaussian 
distribution with mean magnitude —0.26 and a = 0.2. These as- 



sumptions are mostly consistent w ith observation s dStanek et al.l 
and the fitted distribution from lUdalskilfcOOOh . but there may 
be small offset s betwe en local and bulge red clump giants. It was 
noted in Sumi ( 2004) that there is some as-yet unexplained off- 
set (0.3 mag) in the extinction-corrected mean RCG magnitudes in 
the OGLE fields. A possible explanation for this offset is that the 
RCG population effects are large: so that the absolute magnitude 
of RCG stars is significantly differe nt for RCGs in the bulge com- 
pared to local RCGs, as claimed bv lPerciva l & Salaris] J2003b and 
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Figure 5. Cross-correlation term Cib for 45 OGLE-II Galactic bulge fields 
as a function of field Galactic longitude (top) and latitude (bottom). Open 
circles in the top plot of Cib vs. I correspond to fields 6, 7, 14, 15 and 43 
which have relatively extreme galactic latitudes, see Fig.[T] 



ISalaris et alj d2003l) . A different value of the distance to the Galac- 
tic centre to that assumed here (8 kpc) would in part account for 
the discrepancy, however would not remove it completely. Usin g a 
value of 7.6 kpc aEisenhauer et al . 2005; Ni shivama et al.l2006h as 
the distance to the Galactic centre would change the zero-point by 
0.12 mag, resulting in an offset value of 0.18 mag. It is also possi- 
ble that reddeni ng toward th e Galactic centre is more complicated 
than assumed in lSumil 1 2004). In order to compare the model proper 
motion results with the observed data, it was necessary to shift the 
mean model RCG magnitudes to correspond with that observed in 
each of the OGLE fields. The model RCG magnitudes were fitted 
with a Gaussian curve. The mean of the model RCG magnitudes 
was then shifted by a value Am, see Table [3] to correspond with 
the observed mean RCG magnitude in each of the OGLE fields. No- 
tice that we concentrate on second-order moments (proper motion 
dispersions) of the proper motion, so a small shift in the zero-point 
has little effect on our results. 

Every model particle has an associated weight, w%. The par- 
ticle weight can take values < Wi < 20. In order to account 
for this weighting, \wi] stars are generated for each particle with 
the same kinematics but magnitudes determined as above. \wi\ is 
the nearest integer toward +oo. Each model star is then assigned 
a weight, 7* = tUj/puii]. Notice this procedure allows us to in- 
crease the effective number of particles to better sample the lumi- 
nosity function. The total number of stars and the number of stars 
assigned RCG magnitudes in each field are listed in Table[3]as n a ii 




Figure 6. Galactic kinematics from the model of Debattista et al (2007, in 
preparation). Bulk stellar motion in the mid-plane of the Galaxy is shown 
super-imposed on the stellar density. The Sun is located at the origin (not 
shown). An example line of sight is shown. The model can be rotated to four 
equivalent positions for each line of sight due to symmetry (see section |3~TV 



3000- 




Figure 7. Number count as function of apparent magnitude, /, for OGLE-II 
field 1 . The number count histogram is shown along with the fitted function 
equation {6}. The fraction of RCG stars, R^, is evaluated over the magni- 
tude range 7 p ± 3cr for each of the (fc = 1 . . . 49) OGLE-II fields. The ratio 
Rk is assumed to be the same at all stellar distances for each field. 



and n rcg respectively. 81806 stars from the model were used to 
compare model kinematics to observed values. 



3.2 Model kinematics 

Stars with apparent magnitudes within the limits m m i n = 13.7 and 
Ti max = 15.5, were selected from the model data. This magnitude 
range corresponds to the selection criteria imposed on the observed 
data sample, see section |2~2l Model stars with total proper motions 
greater than 10 masyr -1 (corresponding to > 380 kms -1 at a 
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distance of the Galactic centre) were excluded on the basis that 
such stars would be similarly excluded from any observed sample. 
The fraction of weight removed and number of stars removed in 
this way only amounted to a few per cent of the total weight and 
number of stars in each field. Bulge model stars were selected by 
requiring a distance d > 6 kpc. 

The mean proper motion and proper motion dispersions in the 
latitude and longitude directions were computed along with their 
errors for all model stars in each field which obey the above selec- 
tion criteria. The weights on model stars, 74, were used to compute 
these values. 

We then tested whether the finite and discrete nature of the 
model data gives rise to uncertainties in the measured proper mo- 
tion dispersion values. We measured the intrinsic noise in the model 
by comparing the proper motion dispersions computed for four 
equivalent lines-of-sight through the model for each field. The 
spread of the proper motion dispersions for each field was then used 
as the estimate of the intrinsic noise in the model. The mean (me- 
dian) value of these errors in the longitude and latitude directions 
are 0.08 (0.06) and 0.12 (0.097) masyr" 1 respectively. 

The statistical error for the proper motion dispersions in the 
longitude and latitude directions for each field were combined in 
quadrature with the error arising from the finite discrete nature of 
the model data to give the total error on the proper motion disper- 
sions computed from the model. 



4 COMPARISON BETWEEN THEORETICAL MODEL 
AND OBSERVED DATA 

The observed and predicted proper motion dispersions for each of 
the OGLE-II fields are shown in Table|4] Fig.[8]shows the observed 
proper motion dispersions for each of the analysed OGLE-II fields 
plotted against the predicted model proper motion dispersions. 

Fig- 00 shows that the model predictions are in general agree- 
ment with observed proper motion dispersions for the OGLE- 
II fields. The model has been used previously to predict the 
proper motion dispersi ons of 427 starsQ entries observed by 
ISpaenhauer et al. 1 19921) in a single 6' x 6' field toward the bulge 
dBissantz et alj|2004l) . The model value of a\ in this previous anal- 
ysis was in agreement with the observed value, yet the model and 



The significant difference between the observed proper mo- 
tion dispersions of adjacent fields (e.g. fields 1 and 45) might hint at 
some fine-scale population effect, whereby a group of stars surviv- 
ing the selection criteria have a significant and discrepant kinematic 
signature. Higher- accuracy observati ons using the HST supp ort this 
evidence of such population effects dKoztowski et al.ll2006h . 



observe d values of Ob were signifi cantly different. The 6' x 6' field 
used by ISpaenhauer et alj dl992h falls within the OGLE-II field 
number 45. The model prediction of o\ for stars in OGLE field 45 
is completely consistent with the measured value. The model pre- 
di ction of Ob shows a s imilar discrepancy to the previous analysis 
of iBissantz et alj d2004l) . 

Fig. [9] shows the ratio 7? = ai/a^ and cross-correlation term 
Cib = o"ib/(<7i<7b) computed using the model and observed data. 
Typically the model predicts more anisotropic motion with R > 1 
than what is observed. 

The model predictions for stellar kinematics in the latitude di- 
rection may be problematic. This is not surprising as the model is 
not well constrained toward the plane due to a lack of observational 
data because of the heavy dust extinction. The problem is currently 
under investigation. Similarly, the model predictions for ai degrade 
as I increases. This is because the model performance has been op- 
timised for regions close to the Galactic centre. 



There are two repeated entries in Table 2 o flSpaenhauer et alj ll992l) . 



No attempt has been made to account for the blending of flux 
inherent in the OGLE-II crowded-field photometry. It is certain 
that a fraction of st ars in each OGLE-II field suffers from some 
degree of blending dKoztowski et alj|2006h . To investigate this ef- 
fect, we checked one field c overing the lens M ACHO-95-BLG- 
37 (I = 2.54°, b = 3.33°, iThomas et alj|2005h from the HST 
proper motion survey of lKoztowski et al.1 d2006t) . which falls inside 
OGLE-II field number 2. HST images suffer much less blending, 
but the field of view is small, and so it has only a dozen or so clump 
giants. We derive a proper motion of a\ = 3.13 ± 0.57 masyr -1 , 
andob = 2. 17±0. 40 mas yr _ . These values agree with our kine- 
matics in field 2 within 0.2cr for o\ and 1.6a for <7b- The errors in 
our proper motion dispersions are very small (~ km s _1 at a dis- 
tance of the Galactic centre), but it is likely that we underestimate 
the error bars on the observed data due to systematic effects such 
as blending. 



4.1 Understanding the differences 

We now seek to understand the cause of the differences between 
the model and the Milky Way, at least at a qualitative level. We 
first consider the possibility that the difference can be explained by 
some systematic effect. We compute the differences between ob- 
served proper motion dispersions of nearest fields for fields with 
separations less than 0.25 degrees. No pair of fields is used twice, 
and the difference A = a, — aj is always plotted such that 
> \bj\. Ai >0 b s and Ab,obs denote the difference in observed 
proper motion dispersions between adjacent fields in the longitude 
and latitude directions respectively. The equivalent quantities pre- 
dicted from the model are denoted Ai. mo( j and Ab.mod- In Fig. 1101 
we see that the deviations Ai i0 b s and Ab,obs scatter about 0, but 
have a quite broad distribution in both the I and b directions, with 
several fields inconsistent with zero difference at la (defined as the 
sum in quadrature of the uncertainties of the corresponding quanti- 
ties of the two fields under comparison). Several deviations are as 
large as 0.2 mas yr -1 (corresponding to ~ 8 kms -1 at the Galac- 
tic centre) and many a's away from zero. In view of the fact that 
these differences have mean close to zero, it is possible that these 
deviations are due to some systematic effect rather than to intrinsic 
substructure in the Milky Way. We return to this point briefly in the 
discussion. 

In the case of the model uncertainties, however, Fig. II llshows 
that in most cases the differences Ai >mo( j and Ab, mo d are consistent 
with zero at the la level, indicating that these error estimates are 
robust. 

We now seek to explore the correlations of the residuals with 
properties of the model. We plot residuals 5i,b = (cmod — fobs), 
where <r mo d and a Q b s are the model and observed proper motion 
dispersions in the corresponding Galactic co-ordinate. The errorbar 
length is (u^ od + u^) 1 ' 2 where w mo d and u hs are the uncertain- 
ties in the model and observed proper motion dispersions, respec- 
tively. Plotting these quantities as a function of I, we note that there 
is no significant correlation, but that the largest deviations in the 
latitude proper motion dispersion occur close to I — 0, see Fig. 1121 
In plotting <5i,b as a function of b, the reason which becomes evi- 
dent is that the fields closest to the mid-plane have the largest 8b, 
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Model <7i (mas yr ) Model cr b (mas yr ) 

Figure 8. Comparison between observed and predicted proper motion dispersions for stars in the OGLE-II proper motion catalogue o flSumietalJj2004l) . Left: 
Proper motion dispersions in the galactic longitude direction, o\. The OGLE-II field number is indicated adjacent to each point, see also Fig. [T] Fields with 
galactic longitude |/| > 5° are shown in magenta; fields within Baade's window are shown in red; all other fields in blue. Right: Proper motion dispersions in 
the galactic latitude direction, <r b , shown with the same colour scheme. 




Model <Ti/<7b Model Gb 

Figure 9. Left: Ratio of proper motion dispersions R = <T\/a-^ for the observed OGLE-II proper motion data and model predictions. The model generally 
predicts more anisotropic motion, i.e. R > 1 than is observed in the data. Right: The cross-correlation term Cib = (Tib/ciov 



see Fig.Q~3] The density distribution in this region is uncertain due 
to presence of dust and the large extinction corrections required. 
This may explain why the residuals of <7b seem to correlate more 
with b than those of a\. We note that the <t\ residuals also seem to 
have some dependence on 6. A possible explanation is that there is 
some additional effect due to dust which has not been accounted 
for. 



5 DISCUSSION 

Red clump giant stars in the dense fields observed by the OGLE-II 
microlensing survey can be used as tracers of the bulge density and 
motion over a large region toward the Galactic centre. The proper 
motion dispersions of bulge RCG stars in the OGLE-II proper mo- 



tion catalogue of lSumi et alj J2004I) were calculated for 45 OGLE- 
II fields. The kinematics derived from the ground-based OGLE-II 
data were f ound to be in agreeme nt with HST observations in two 
fields from Kozlowski et al.1 bfJOot) . It is reassuring that the results 
presented here are consistent with those derived from the higher 
resolution HST data, despite possible selection effects and blend- 
ing. 

The observed values of a\ and <r b were compared to predic- 
tions from the made-to-measure stellar-dynamical model of De- 
battista et al. (2007, in preparation). In general, the model gives 
predictions qualitatively similar to observed values of o\ and <7b 
for fields close to the Galactic centre. The model is in agree- 
m ent with observed OG LE-II data in the direction previously tested 
bvlBissantz et al.l d2004T) . Using the definition of |Pe Lorenzi et all 
d2007h . the effective number of particles in the model used here 
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Table 4. Proper motion dispersions in the longitude and latitude directions, <ti, <r b , and cross-correlation term Cih for bul ge stars in 45 OGLE-II fields. High 
precision proper motion data for bulge stars were extracted from the OGLE-II proper motion catalogue |Sumi et al. 2004!). N is the number of stars selected 
from each field. Field 44 was not used due to the low number of RCGs in this field. 
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0.03 ±0.08 


-0.12 ±0.01 


9799 


17 


5.28 


-3.45 


3.09 ±0.09 


2.81 ±0.02 


1 ^fl _I_A AT 

Z.JU ±U.U/ 


2.42 ±U.U1 


-0.01 ±0.11 


-0.12 ±0.01 


10268 


18 


3.97 


-3.14 


3.20 ±0.09 


2.92 ±0.02 


1 A O 1 A AO 

2.48 ±0.08 


o ai i r\ m 

2.62 ±0.02 


0.02 ±0.03 


-0.13 ±0.01 


14019 


19 


4.08 


-3.35 


3.06 ±0.13 


2.90 ±0.02 


1 A A _l_ A 1 6. 

ZAy ±U.zo 


2.6U ±U.U2 


0.01 ±0.03 


-0.17 ±0.01 


13256 


20 


1.68 


-2.47 


3.12 ±0.06 


3.27 ±0.01 


1 ZZ _I_A AC 

2.66 ±0.05 


o on 1 n n i 

2.82 ±0.01 


0.07 ±0.03 


-0.12 ±0.01 


17678 


21 


1.80 


-2.66 


3.12 ±0.06 


3.31 ±0.02 


1 C*7 _I_A AO 

Z.J I ±U.Uo 


2.yU ±U.U2 


-0.02 ±0.03 


-0.13 ±0.01 


17577 


22 


-0.26 


-2.95 


3.17 ±0.04 


3.17 ±0.02 


1 A £L 1 A 1 1 

2.46 ±0.12 


o A in no 

2.84 ±0.02 


0.01 ±0.03 


-0.01 ±0.01 


19787 


23 


-0.50 


-3.36 


3.13 ±0.17 


3.15 ±0.01 


1 Z1 in 1 A 
Z.OZ ±U. 1(J 


1 o/f ^_n m 
2.84 ±0.02 


-0.02 ±0.14 


-0.04 ±0.01 


17996 


24 


-2.44 


-3.36 


2.77 ±0.04 


2.96 ±0.01 


1 11 I A 1 A 

2.32 ±0.10 


o ^ o n n I 

2.48 ±0.01 


-0.04 ±0.04 


0.02 ±0.01 


16397 


25 


-2.32 


-3.56 


2.76 ±0.07 


2.91 ±0.01 


1 /1 1 1 A I C 

z.4/ ±U.1j 


2.50 ±0.01 


-0.04 ±0.03 


0.02 ±0.01 


16386 


26 


-4.90 


-3.37 


2.80 ±0.17 


2.68 ±0.02 


111 1 A A A 

Z.ZZ ±0.04 


o i -7 n n i 

2.17 ±0.01 


-0.00 ±0.03 


0.02 ±0.01 


13099 


27 


a ai 


a A^ 
-3.65 


2.78 ±0.07 


2.63 ±0.02 


1 1 A _l_ A ("1 /I 

ZAy ±U.U4 


2.15 ±0.01 


-0.06 ±0.02 


0.03 ±0.01 


12728 


28 


-O. /O 


-4.42 


3.02 ±0.11 


2.63 ±0.03 


1/1/1 _I_A 1Z 

z.44 ±U.3o 


t n ^n no 
2. 12 ±0.02 


0.05 ±0.06 


-0.01 ±0.01 


8367 


29 


-6.64 


-4.62 


3.02 ±0.21 


2.66 ±0.03 


1.79 ±0.14 


2.09 ±0.02 


-0.00 ±0.11 


-0.02 ±0.01 


8108 


30 


1.94 


-2.84 


3.13 ±0.07 


3.04 ±0.02 


2.59 ±0.11 


2.70 ±0.02 


-0.04 ±0.08 


-0.12 ±0.01 


17774 


31 


2.23 


-2.94 


3.08 ±0.05 


3.11 ±0.02 


2.68 ±0.11 


2.74 ±0.01 


0.08 ±0.05 


-0.12 ±0.01 


17273 


32 


2.34 


-3.14 


3.10 ±0.09 


3.10 ±0.02 


2.56 ±0.04 


2.78 ±0.01 


0.11 ±0.02 


-0.13 ±0.01 


15966 


33 


2.35 


-3.66 


2.82 ±0.11 


3.08 ±0.02 


2.57 ±0.12 


2.77 ±0.02 


0.07 ±0.06 


-0.14 ±0.01 


15450 


34 


1.35 


-2.40 


3.18 ±0.06 


3.36 ±0.02 


2.62 ±0.03 


2.92 ±0.01 


0.04 ±0.02 


-0.11 ±0.01 


16889 


35 


3.05 


-3.00 


3.05 ±0.05 


3.09 ±0.02 


2.59 ±0.07 


2.72 ±0.02 


0.08 ±0.03 


-0.14 ±0.01 


15973 


36 


3.16 


-3.20 


3.00 ±0.06 


3.19 ±0.02 


2.95 ±0.40 


2.77 ±0.02 


-0.05 ±0.08 


-0.16 ±0.01 


14955 


37 


0.00 


-1.74 


3.29 ±0.04 


3.29 ±0.02 


2.70 ±0.04 


3.04 ±0.01 


-0.01 ±0.01 


-0.05 ±0.01 


20233 


38 


0.97 


-3.42 


3.01 ±0.07 


3.15 ±0.01 


2.60 ±0.14 


2.84 ±0.02 


0.07 ±0.07 


-0.12 ±0.01 


15542 


39 


0.53 


-2.21 


3.22 ±0.03 


3.21 ±0.01 


2.69 ±0.06 


3.00 ±0.01 


0.01 ±0.04 


-0.07 ±0.01 


24820 


40 


-2.99 


-3.14 


2.84 ±0.04 


2.84 ±0.01 


2.28 ±0.07 


2.47 ±0.02 


-0.09 ±0.07 


0.05 ±0.01 


13581 


41 


-2.78 


-3.27 


2.86 ±0.06 


2.78 ±0.01 


2.60 ±0.19 


2.41 ±0.02 


-0.16 ±0.07 


0.04 ±0.01 


14070 


42 


4.48 


-3.38 


3.07 ±0.05 


2.89 ±0.02 


2.44 ±0.15 


2.63 ±0.02 


0.02 ±0.02 


-0.15 ±0.01 


10099 


43 


0.37 


2.95 


3.13 ±0.06 


3.17 ±0.02 


2.72 ±0.10 


2.87 ±0.01 


0.04 ±0.07 


0.02 ±0.01 


11467 


45 


0.98 


-3.94 


3.02 ±0.05 


2.97 ±0.04 


2.42 ±0.14 


2.61 ±0.04 


0.06 ±0.11 


-0.13 ±0.02 


2380 


46 


1.09 


-4.14 


2.87 ±0.08 


2.90 ±0.04 


2.53 ±0.21 


2.67 ±0.04 


-0.03 ±0.06 


-0.16 ±0.03 


1803 



is 3986. This relatively low number results in large errors on the 
model proper motion dispersions and we therefore recommend re- 
garding interpretations based on this model with some caution. An 
i mproved model with a lar ger number of particles (the recent study 
bv (Pe LorenzTe t al. (2007) has an effective particle number ~ 10 6 ) 
will decrease the errors on the model predictions and allow a more 
useful comparison between model and observed proper motion dis- 
persions. 

The OGLE-II fields mostly extend over ~ 17° in longi- 
tude and about 5° in latitude across the Galactic bulge region 
and can therefore provide a more powerful set of constraints on 



stellar motions predicted by galactic models. The high-accuracy 
proper motion data for the 45 fields and those obtained with HST 
dKozlowski et alJ 2006) can be used as direct input in the made-to- 
measure method to construct a better constrained dynamical model 
of the Milky Way. 

The statistical errors of our proper motion dispersions are 
small (~ kms" 1 ), but systematic uncertainties (for example due 
to incorrect dust extinction treatment) which were not included in 
the analysis may be significant. Nevertheless, it is interesting to 
note that there appears to be significant difference between the ob- 
served proper motion dispersions of adjacent fields (e.g. fields 1 
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Figure 10. Difference between observed proper motion dispersions for 
pairs of fields with separations less than 0.25 degrees (corresponding to 
~ 40 pc at the Galactic centre). 

and 45). This might hint at some fine-scale population effect, where 
the kinematics of the bulge may be not in total equilibrium (e.g. due 
to a small accretion event). Higher-accuracy observations using the 
HST ma y provide further e vidence of such population effects. We 
note that lRich et al.l {2006) suggest the possible existence of cold 
structures using data from a radial velocity survey of Galactic bulge 
M giant stars although their conclusion could be strengthened by a 

larger sample of stars. 

The OGLE-II proper motion catalogue JSumi et alj[2004h for 
millions of bulge stars is still somewhat under-explored. For exam- 
ple, it will be interesting to explore the nature of the high proper 
motion stars (p > 10 masyr -1 ) and search for wide binaries in 
the catalogue. Some exploration along these lines is under way. 
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